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How can we tell whether dark energy is composed by multiple fields? 
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Dark energy is often assumed to be composed by a single scalar field. The background cosmic 
expansion is not sufficient to determine whether this is true or not. We study multi-field scalar-tensor 
models with a general dark matter source and write the observable modified gravity parameters 
(effective gravitational constant and anisotropic stress) in the form of a ratio of polynomials in the 
Fourier wavenumber k of order 2N, where N is the number of scalar fields. By comparing these 
observables to real data it is in principle possible to determine the number of dark energy scalar 
fields coupled to gravity. We also show that there are no realistic non-trivial cases in which the 
order of the polynomials is reduced. 


I. INTRODUCTION 


Testing the nature of dark energy and its possible modification of gravity at very large scales is currently one of the 
most interesting research activities in cosmology. Modifications of standard gravity are often modeled by introducing 
additional mediating fields in the gravitational Lagrangian. One of the most well studied example is the so-called 
Horndeski theory, which adds to the Einstein-Hilbert Lagrangian a single scalar field that obeys the most general 
second order equation of motion [l[. 

It is however in principle possible that dark energy is actually a composite state of two or more scalar fields 
interacting through a common potential and coupled to gravity with different strengths (see e.g. HI). Since the 
only observable at background level is the Hubble function H{z) as a function of redshift, and since to any H{z) one 
can always associate a particular single-field potential (see appendix), it is impossible to tell whether dark energy 
is driven by multiple fields instead of just one, as almost universally assumed. In this paper we show how linear 
cosmological perturbations produce instead distinguishable effects that depend on the number of fields. 

As shown in several papers (e.g. [MU) , a generic modification of gravity introduces at linear perturbation level two 
new functions that depend only on background time-dependent quantities and, in Fourier space, on the wavenumber k. 
One function, that we denote here with Y{t, k) (sometimes also called Ges), modifies the standard Poisson equation, 
while the second one, rj{t,k), the anisotropic stress or tilt, provides the relation between the two gravity potentials 
dt, $. In standard gravity, one has Y = ry = 1. 

In the so-called quasi-static regime (i.e. for linear scales that are below the sound horizon) the two functions V, ? ? 
take a particularly simple form in the Horndeski models and can be directly constrained through observations [ij, [l^ . 
This also holds true in some cases (see also [1^ and [I^) of bimetric models [13. 

Since the equation of motion of the Horndeski Lagrangian is second order, the equations contain at most second 
order space derivatives and therefore terms in Fourier space. It is then no surprise that the resulting forms of the 
functions Y, rj include indeed polynomials of second order in k 


fi + ehi\ _ /i + phsX 

^ dt ^yi + k'^h^) ’ AirGa^pS ^\l + k‘^h^) 


( 1 ) 


where the hi functions are time dependent functions that depend on the specific Horndeski model, a is the scale factor, 
p is the average matter density in the Universe, 5 is the matter density contrast and $, are the two perturbation 
potentials, to be defined below . Also the combination Y{\ + p) that appears in the lensing potential equation has 
the same structure. As shown in e.g. [l^, the fact that the k structure is fixed regardless of the specific model 
(within the Horndeski or the bimetric gravity class), allows one to combine observations of weak lensing, redshift 
distortions and galaxy clustering to constrain or detect modifications of gravity at cosmological scales in a relatively 
model-independent way. It is therefore interesting to inquire about how general the form of Eq. CD is. 

In 0 (see also [ll|) it has been shown that models that obey n-th order equations of motion generalize the Re¬ 
structure of Y, p so that 


. Pn\k) ,, , P^^\k) 

V = ^ 2 — 7 ^ 77 — , y = ^ 1 - 


Pi^\k) 


Pi'^Hk) 


( 2 ) 
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where Pn\P^\P^'^ are even polynomials of n—th order in k with time-dependent coefficients, normalized for 
convenience so that Pn^\o) = P^\o) = Pn^\o) = 1. This case includes Lagrangians which depend on several scalar 
fields, because they also obey a system of second-order differential equations that is dynamically equivalent to a single 
higher-order equation. A Lagrangian built out of several scalar fields is in general stable provided the interaction 
potential is bounded from below. 

In order to compare the models to observations however one needs the explicit form of the polynomials: this is 
the main goal of this paper. Here we consider the simplest non-trivial scalar tensor theory, namely a multi-field 
Brans-Dicke model, parametrized by several coupling constants uji and by a general potential (/) 2 ,...). Moreover, 

we generalize the Brans-Dicke model also in a different way: we include a generic fluid matter source with equation 
of state w and sound speed c^, both of which generally time dependent. This might turn out useful to compare real 
data to models in which dark matter includes hot or warm components. 

Finally, we investigate whether there are cases in which the higher-order k terms in the polynomials in Eq.([2]) cancel 
out, thereby effectively reducing the order of the system and, consequently, the number of observable scalar fields. We 
find that this is not possible unless the fields are decoupled from gravity or the Universe is filled only with radiation. 
We argue therefore that the fc-structure of Y,r] is a unique signature of the number of dynamically active dark energy 
fields coupled to gravity. 

Multi-scalar-field scenarios are being considered in gravitational physics and cosmology since a while. In the context 
of modified theories of g ravity such models are studied in e.g. [I9l - l22ll. for string theory inspired multifield dark energy 
models see e.g. [H, [ 2 J while for the inflationary context see e.g. |25l - [2^ . 


II. STANDARD BRANS-DICKE THEORY 

Let us first derive the cosmological perturbation equations for the standard Brans-Dicke theory with a potential. 
The corresponding action is: 


5'i = J - V{cj))) + j (fxy/^C matter i (3) 

where as matter fluid we consider a general component characterized by an equation of state w and a sound speed Cg, 
both arbitrarily time dependent. The equation of motion (e.o.m.) emerging from the variation with respect to the 
metric is 


+ [n</' + 


+ ^5;.. = SttGT^ 


(4) 


where 


= (p + p)u^Uu + pg^v (5) 

is the stress-energy tensor of the perfect fluid, with p being the energy density, p being the pressure and the 
4-velocity. We will need the trace of Eq. Q which reads as 

-(t)R + 30(1) + 2^ (V<^)2 - -b 2V{(I>) = SttGT. (6) 

(j) (j) 


The scalar field e.o.m. is 


R + 2 


^BD 

(t> 


□ (/) 


^BD 



- U0 = 0, 


where the derivative with respect to the field (j) is denoted as U^. Using Eq.® we obtain 


(7) 


{3 + 2ujBD)a(t> + 2V{(l))-V^ = 8nGT. (8) 

Erom this equation it follows that in the limit ujbd 00 the scalar field is constant (see e.g. [ 2 ^) and therefore 
in this large-uiBD limit the gravitational equation of motion in Eq.® coincides with the Einstein-Hilbert equation of 
motion. It is interesting to note that the Brans-Dicke Lagrangian reproduces also the so-called f{R) models in the 
limit in which ujbd vanishes and 4> = df{R)/dR (see e.g. [30|). 
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We start deriving now the cosmological perturbation equations for the Brans-Dicke theory in a spatially flat 
Friedmann-Lemaitre-Robertson-Walker (FLRW) metric. We assume a single fluid matter source with general equation 
of state and general sound speed. As it is well know, at first order in the perturbation parameter one can bring the 
scalar line element into the form 


= a^(r)[—(1 -|- + (1 + 2^)dx’'dxi], (9) 

where r is the conformal time, a is the scale factor and 'k, $ are the scalar potentials. We also decompose the scalar 
field into the background sector and the perturbed sector as x) = + (f{t, x). After switching to Fourier space 

we obtain from the (0,0) element of Eq.(|4|) 

= ( 10 ) 

(p 1 

where 'H. = a! ja is the conformal Hubble function, the prime denotes differentiation with respect to the conformal 
time r and dp denotes the matter density perturbation. 

Now we assume the validity of the so called quasistatic approximation (see e.g. [l^ IMI), which means we con¬ 
sider only subhorizon scales 'Pi} and also that the perturbation fields vary slowly enough with time so that 

$, 4', 4'', (/?, p' are all negligible with respect to dp unless multiplied by or by the effective scalar field mass 
NP = H 00 . Note that in case of large ojbDi the terms like the last term on the first line of Eq. (flUl) is not very large 
because all such terms multiply the derivative of the field, which is always of order of XjuiBOi as we have discussed 
above. In this approximation the previous expression becomes 

2(j)—^-\—-(f — SnGdp = 0. (11) 


Perturbing the scalar field e.o.m. we obtain 


- 12 {n' + ■H^)^'] -k - 6 ^(^'' - 3 $') - + 

, r , a ; BD „ 0 v 7 2 ujbd ,,2 , 2 ujbd , ,, ,, „ 

2——V Vip -k 2——V''Vo<k>- = 0. 

(p (p P P 4> 


In the quasistatic limit this becomes: 


( 12 ) 


.P ^ ^ ,^UJBD k , 2 


4—$-k2—^'-(2- 


-kM^)^ = 0 . 


' (j) P 

Finally we use the traceless part of the gravitational e.o.m. to derive one more perturbation equation: 

(j>Gl - 3V°Vo<^ - 2V^V,4> + 4>R- - ^v{,4>) = 87rG(r° - T). 

2 (f) 2 (f 2 

Recalling that dTp — dT = —dp — (3c^ — l)dp = —3cldp , where = dp/dp is the sound speed, we get: 

nj2 nj nj2 9 ‘X h'^ 7-/2 'IJ' 

- + <^[- 64 $' + e^vi/ _ + 4<k>" + 2^(k’ + 

-\(j)[-2kH - -k 12('H' -k n^)m + - 3$') 

+ - 3^VVVo<(> - = -2A^Gcldp. 

2 (b^ 0 2 


In the quasi-static limit this reduces to 


7 2 A; ^ 

2(^—4) -k 2—p -k 2cj)—A! = —2A'!TGcldp. 


(13) 


(14) 


(15) 


( 16 ) 
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Collecting all the above-derived perturbation equations together we have the following system of equations: 


7 2 

2(/)—$ -I- 2—if + 2(/)—vh = —24:TrGc^Sp, 


4 $ + 2t!»-(2i^2£^l + M> = 0, 


^2 jpi 

2(j)—^ H— -5ip = SirGSp. 


(17) 

(18) 
(19) 


In case of = 0, these expressions agree with Q. Solving the system and replacing Sp with 6 = Sp/p = 
8ttG{Sp )we obtain: 


= - 


3VLm'H‘^8 




( 20 ) 


where the modified gravity parameter Y is defined to be 


1 (1 -|- 3c^)M^(j} + fT(4 -l- 2ujbd + 64(1 -I- cobd)) 
<!> -f ^(3 -|- 2ujbd) 


Let us calculate the first two correction terms on top of the GR case. Since local gravity measurement constrain the 
Brans-Dicke coupling parameter to be very large , the expansion of Y up to second order in powers of 1 /wbd gives 


Y - 


(1 + 34) 


1-34 

2ujBD(p 


(34 




( 22 ) 


Note that in order for the coupling constant in the zeroth order Poisson equation to be equal to Newton’s constant 
at the present time, we should require (()(0) = 1. We define the anisotropic stress as (see e.g. El) 


< j ) -j- ^(2 -|- 3 Cg -|- 2lobd) 

(1 + 3cf)M^p -\- ^(4 - l - 2ujbd + 6 cg(l - I - ujbd)) 


If we set = ujbd = 4 = 6 , we get + = 4/3(^, p = 1/2, which, as expected, are the same values as in a standard 
f{R) model in the limit of small mass M provided that (p = df{R)/dR. By comparing those expressions against 
observations one can in principle estimate the coefficients in the polynomials and constrain the values of the matter 
sound speed and the field mass. However, we can give a rough estimate for the numerical values of the coefficients 
at the present epoch; for that purpose we use that ^( 0 ) = 1 from the above discussion, assume the potential to be 
of the form V{(P)/16ttG = /32'kG and consider the case of 4 0. From the Friedmann equation we then have 

^ V{p) and therefore ^ 6774 where Hq is the present value of the Hubble function. Using those and the 
fact that at the present epoch we have a = 1, 77 and Y can be written as: 

(24) 


(25) 

The second order expansion of p in IjujBD is 


Y = 


1 -f 


p = 


2k^ 


1 + 
1 + 


fc 2 




1 + 


2k^ 


1 


p ~ 


1 + 34 


(34-1)(34 + 2) 
2(34 + 1)"wbb 


+ ( 9 Cs -b 34 


2 ) 


(1 -b 3c;)M^p -b (4 -b 64) 


.2\k-‘ 


4(34 + 1 ) 


3 2 

TtUJ 


BD 


(26) 


Notice that from Eqs. (EU and it follows that for cf = 1/3, i.e. for a relativistic fluid, both Y and p are 
independent of ujbd- This reflects the fact that under a conformal rescaling the Brans-Dicke theory can be recast 
into a ordinary gravity model with field-matter coupling proportional to the trace of the energy momentum; when 
this vanishes, as for a relativistic fluid, the coupling vanishes as well. 
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III. TWO-FIELD BRANS-DICKE THEORY 

We move now to the case in which the Lagrangian includes two fields, both coupled to gravity. We adopt then the 
following two-field Brans-Dicke Action 


S 2 = j d^xy/^[(l)iR- d^(j)id^(t)i +(I>2R- ^g^''d^(j)2du<t)2 - (^ 2 )] 

-I- J d'^x^/^C 

matter ■ 

The gravitational e.o.m. of this theory is 


(27) 


(^1 + 4>2)Gij_v + [□((/'I + 4>2) + 2 

-V^V,(<(.1 + d)2) - = SttGT^,. 

(pi (p2 I 


Taking the trace we obtain 


- (<(>1 + 4>2)R + 3n(<^i + (t>2) + + '^iy<t>2f + 2V{(i>i,h) = SttGt. 

01 02 


The 01 and 02 e.o.m. are 


R + 20n<?5>i - ’^d°'(i)ida(j)i - = 0, i = 1,2. 


(28) 


(29) 


(30) 


Summing Eos. d^ . dHni) together we obtain: 

(3 -l- 2uji)0(j)i -|- (3 -|- 2u!2)Gi(j)2 + 2Y(())i, (f>2) — ~ ^,<1)2 ~ SttGT. 


(31) 


Repeating the arguments after Eq.® we conclude that this theory reduces to standard general relativity in the 
large-wi,W 2 limit. 

Proceeding as in the previous sectiom, the perturbed (0,0) component of Eg. (1281) in the quasi-static approximation 
leads to the corresponding perturbation equation: 


^2 ^2 

2(^1 + 02)^d) -I- -it{(Pi + 1 P 2 ) - SnGSp = 0. 

Similarly, from the scalar field e.o.m. Eas. (l30l) we derive the corresponding perturbation equations: 


(32) 


4—$ + 2—[2^—=0, i = l,2 

0i 


(33) 


where Mf = i = 1,2. Perturbing also the (0,0) component of the traceless part of the gravitational e.o.m. in 

the quasistatic limit, we obtain the following system of perturbation equations: 


2((/'i -I- ((' 2 )^$ + + P 2 ) + 2(</)i -I- 4>2)^'^ = -2AttGc16p, 


2 vr J- ' ' "Vt-L ' 2 

U2 . 

t 2 /-,2 /A./i 2 


4^$-h2^4' - [2^^ + M'^]pi = 0,i = 1,2 

(pi 

^2 ^2 

2[(j)i -I- (j) 2 )^^ + + P 2 ) = SnGSp. 

Solving this system and replacing Sp with S = Sp/p = 8'!rG{6p)/3flmH'^ as before, we obtain finally 


(34) 

(35) 

(36) 




2 


( 37 ) 





where 
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. 1 + A2k^ + A3k^ 

“ h + y44fc2 + ^5^4 ’ 


with the coefficients being 


Ai 

A 2 


^3 

A 4 

A 5 


l + 3c4 
(/>! + 02 ’ 

2 3Cg)(0i + 02)wi + (2 + 3Cg) 01 

(3c2 + 1) (01 + 02)a2 Mf0i 

(1 + 3c^) (01 + 02)a;2 + (2 + 3Cs) 02 

M|02 

4(2 + 3Cg)(w20i + wi02) + 4wia;2(0i + 02)(1 + Sc^) 

(1 + 3c2) M^M|0 i 02(01 + 02)a4 

2a;i(0i + 02) + 301 2 w2(0i + 02) + 302 


(01 + 02)a2 


M2<(,1 


M202 


4a;ia;2(0i + 02) + 6(0iW2 + 02Wi) 
M2M|0i 02(01 +02)a4 ■ 


(38) 

(39) 

(40) 

(41) 

(42) 

(43) 


Considering the limit of uii and W 2 both much larger than unity and keeping the first terms of expansion which are 
non-trivial in we find for the effective gravitational constant the result: 


Y 3Cs + 1 (1 — 3Cs)(a;i02 + O-’20l) ,^^2 _ (01 + 02) (-^l‘^20j + -^1^102) + 3^(0:102 + ^20l)^ 

01 +02 2aJiW2(01 + 02)^ * 4^w2u;|(0i + 02)3 

Similarly, the anisotropic stress is 


(44) 


^ ^ 1 + g2fc" + Bsk^ 

~ ^^1 + Bik^ + B^k* ’ 

with the coefficients defined as 
' (3c2 + l)’ 

_ 1 r ((3Cs + 2 ) 02 + 2w2(01 + 02)) (01 (3Cg + 2) + 2wi(0i + 02)) 

^ (01 + 02)0^ M|02 Mf 4 >i 

_ 2 (2 + 3c^) (u;i02 + a;20i) + 4a;ia;2(0i + 02) 

^ M|0102(01 + 02)a4 

n _ . __2_ (1 + 3Cg)(0i + 02)a;i + (2 + 3Cg) 0i 

^ ^ (3c2 +1) (01 + 02)a2 

(l + 3Cs) (01 + 02 ) 0:2 + (2 + 3Cg) 02 
M|02 

_ _ 4(2 + 3c^)(o201 + 0102) + 4oiO2(01 + 02)(1 + 3c^) 

® ^ (1 + 3c3) MfM|0i02(01 + 02)a4 

As expected, we find again that + and rj are independent of 01^2 if = 1/3. Notice also that since A 2 

A 3 = B 3 , the lensing potential ih — $ obeys the equation 


(45) 

(46) 

(47) 

(48) 

(49) 

(50) 
B 4 and 


in which the combination +(1 + 77 ) has in general the same fc-structure as Y and 77 . However, for the particular case 
here studied, i.e. the simple Brans-Dicke model, the null geodesics must remain the same as in General Relativity 
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and therefore Y(\ + r]) should be independent of k. In fact we find y(l + 77 ) = (2 + 3c^)/((^i + (/) 2 )- This shows that 
from weak lensing alone it would be impossible to detect the presence of several Brans-Dicke fields. 

In the large uJi and UJ 2 limit, again keeping the first terms which are non-trivial in k'^/a'^, we obtain for rj 


^ 1 (3c^ - l)(3c^ + 2)(a;i^2 + ^ 2 (/>i) 

3c2 + 1 2 (3c2 + 1)^ wia;2(</i'i + ^ 2 ) 

fq 4 I 3 2 _ 2 ) (3Cs + 1 ) + 4>2) + 2 (3Cg + 2) fT(^i</>2 + . 

^ ’ 4(3cHl)3f^a;?a;|(</>i+<^2)2 

It is also interesting to compare our two-field results with the standard single-field ones in the limit when one of the 
Brans-Dicke parameters is very large, i.e. when the corresponding field effectively decouples from gravity. Keeping 
only the zeroth order term in such an expansion and setting (j )2 to zero we have: 


Y - 

r] 


(1 -|- 3c^) 2^ (3Cg(l -|- uji) -|- wi -|- 2^ 

(/)i + ^(2uji + 3)) 

-b ^ (2 -|- 3 Cg + 2a;i) 

(1 -b 3Cg) -b 2^ (3Cg(wi -b 1) -b Wi -b 2) 


(53) 

(54) 


As we see they identically coincide with Eas. d^T]) and In this limit, therefore, the two-field Lagrangian effectively 
reduces to a single-field one. 


IV. MULTI-FIELD BRANS-DICKE THEORY 

Finally, we can generalize the previous steps to the V-field Brans-Dicke. By analyzing the previous perturbation 
equations we can immediately write down the V-field quasi-static perturbation equations as follows: 


^ Z^2 u2 ^ ^ u2 

2 + 2 ^ = -2AttGcISp, 

i—1 i—1 

4^$ + 2 ^ 4 / - [2^^ + = 0, 7 = 1, 2,..., V, 

(pi 

^ ^2 p ^ 

2 ^ 4>i^^ + ^ Z = SirGSp. 

i—1 i—1 


(55) 

(56) 

(57) 


From Fqs. (ES) we have 


Af ,2 

Z^*=2^(2<J> + 'I')Z^^ 

i=i “ 


M? 


(58) 


Therefore Fas. ((551) and ((571) read respectively 


N 


N 


1 


2 Z + 4^(2$ -b 'I')y^ -r 2 - 

z=l 2=1 (()i a-^ ' I 


N 




-b2Z'/>* —4' = -247rGc^dp, 


2=1 


N 


N 


2 Z <^ 4 *^+ 2 ^( 2 $+ 4 /)^ 


^_/c 2 

2=1 ^ <pi ^ 


■M? 


= 8nGSp. 


(59) 

(60) 










Solving this system one finds 




V = 


Y = 




-\-Mf 


(1 + 3cD Ef=i + 2 (2 + 3c2) ^ Ef=i ^ 


(l + 3c?) Ef=i<^. + 2(2 + 3c2)^Et 


fe2 ^AT 

' 2^5+^" 

N 1 


1 




1 fc2 


The expression for r; can now be written as 


by explicitly specifying the coefficients 


„ ^ 1 + 


(61) 


(62) 


(63) 


Co = 


1 


i + 3cr 


(64) 


and 


N 




^ii ^12 ^id 


n "1 




N 


( 2 +3c?) 2 - 5 : ^ n 


Mf 


2=1 

j^l,...,d-l 

n>22>--->2d-l 

N 


'ri 2 Yid-i 


j^i,ii,i2,---,'id-i 


^’■<=^1(>+=!4)(E«2'‘ E n «: 


AT 


Wi, W,;, W; 


'll <(>12 ‘(’id 


j#il,i2,...,id 


2(2 + 3c?)2-^ E n 


M? 


i=l ij^i. *’■ l/i,il,i2,---,id-l 

i=l,...,d-l 
il>i2>->id-i 


(65) 


( 66 ) 


with d = 1, 2,..., N. Here and in the following expressions by the second summation sign in the second line of Ea. dHbll 
we mean multiple sum over the indices ij, j = 1, ...,d — 1, where each of ij runs from 1 to TV omitting the value i. In 
the case of d = 1 it is meant that there is no summation at all. The products of M? are over index j running from 1 
to N omitting the corresponding indices mentioned under the product signs. C* and D* are given by: 

N N 


N N 

D. = {l + 3cl) (E <(>011^1 ■ 

i=i j=i 


(67) 

( 68 ) 


Similarly, we find 


1+eEap* 


'„ 2 i 


Y = C< 


(69) 


where 
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and 


Co = 


1 + 3cs 


(70) 


1 ^ 

<5<i = ^Kl + 3o?)(E«2‘' E 


^ii ^i 2 ^id 


n "i'+ 


il>i2>->id 




N 


2{2 + 3cI)2‘^-^Y. E 


^ii ^i 2 ^id-l 


i—1 ijz^i^ 

j=l,...,d-l 

ii>i2>--->id~i 

N 


n 


M/l 


= E 


^ii ^i 2 ^id 


n "i + 


N 


N 


il>i2>->id 

Wi, Wi, Wi, 




3(e^02-e e n 


Mf 


(71) 


(72) 


i=l i=l ^ 3d^i,il,i2,---,id-l 

j=l,...,d-l 
n>J 2 >--->J<i-l 

with d = 1,2^..., N. C* and U* are given by: 

Af Af 

a = (l + 3c2) {J2cl3,)l[Mf, 

i=l j=l 

N N 

i=i j=i 

Note that Dd = Cd identically for Vd (see also i). 

We do not have any constraint on the values of the coupling constants, the held masses or the Held trajectories, 
therefore the above coefficients will depend in general on the specihc model and should be determined by comparison 
with the experimental data, which is beyond the scope of this paper. 


(73) 

(74) 


V. REDUCING THE POLYNOMIALS ORDER? 

In this section we discuss under which circumstances the polynomials in 77 and Y effectively reduce to lower order 
ones at all times. If this were possible then the observations (at least at Hrst order in the perturbations) would not 
be able to establish the number of active coupled dark energy Helds, not even in principle. 

From Eas. dCT) and (15^ one can see that the order of the polynomials can be effectively reduced only in one of the 
following cases: a) = 1/3, 6 ) at least one of the coupling constants is inHnitely large, c) one or more Mf —>■ 00 , d) 
one or more (j)i —>■ 0, e) there are i and j such that ujij —>■ 0, f)) there are i and j such that ^ and Mf = aM^, 

where a is a non-zero real number. Note that a = — I is conceptually different from all the other values of a; in 
the former case two summands for each pair i,j disappear from the sum, whereas in the latter case the number of 
summands reduces only by one for each pair i,j. We will comment more on this difference below. 

As we have already discussed before, in case of = 1/3, i.e. radiation, our theory is conformally equivalent to the 
standard general relativity and the Helds (regardless of their number) do not have an observable impact except for a 
time-dependence of the effective gravitational constant. 

The case b) is the limit when the Held with inHnite (or very large) coi is completely decoupled from the Brans-Dicke 
theory (as an example of this reduction see the end of Section mil)- The case c) implies that the Held with inHnite 
(very large) mass is static and just renormalizes the gravitational constant. Case d) is unphysical because it leads to 
a singular Lagrangian. In the case e) the kinetic terms of the i and j Helds are absent from the Lagrangian and the 
(jdiR + (j>jR term can be interpreted as a non-minimal coupling of a single Held: il’ijR = + 4'j)R- 

Case f) is non trivial and we need to check the consistency of this case with the Held equations of motion. To do 
so we substitute the condition /) in the (pi equation of motion (1501) : 
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+ 2a^V°Vo0i + 2a- 


LUi 


'i ^,4>i\(pi—0Jicl>j/aujj 0 : 

Yj 

where is the background Ricci scalar. In order for the above equation to be consistent with the (j)j e.o.m: 


Rio) + 2^V°Vo<^, + 


/o 


= 0 , 


(75) 


(76) 


we need that (pj = acpi (leading to uji = u}j) and This means that depending on whether 

a = —lora^^—1, either these two fields decouple completely from the Lagrangian or we can collect the terms 
corresponding to 4 >i and (pj in a way that they effectively act as one field. 

As an example we can consider a potential without (pi, (pj interaction; 


D 


vi(Pi,(P2, ..., M = + v{(Pi,...,kM, 


k^l 


(77) 


where the coefficients Afc are some constant numbers, H is a natural number and hat represents exclusion. For the 
derivatives of this potential we have: 


l<7ii=0j/a 
^,4>j \4’i=4‘j/a 


D 




i(l>3 


\k-l 


vfc-1 






(78) 

(79) 


(i) L — ^ (i) 

In order for the two derivatives to be equal one needs A), = a A), . It is easy to see that this condition is consistent 
with the proportionality of the field masses: \</>i=^j/a- 

To exemplify the a = — 1 case we study the conditions for A 3 , As and B 3 given in Eas. iHTl H51 H51) to be zero. It 
is easy to see that they are consistently zero in one of the following cases: = ^, wi = ^2 = 0 , at least one of 

and M| is infinitely large, at least one of (pi and (p 2 is zero, or, finally, (pi = —(p 2 and oji = W 2 - In the latter 
case, which can be associated with the a = — 1 situation above, both fields completely decouple from the Lagrangian. 
This sheds more light on the case /) above; the cancellation of higher order terms in this case is possible if a pair of 
helds decouples completely from the Lagrangian (a = —1), or if they are still active but in Lagrangian one can group 
together the terms corresponding to the considered pair of fields in such a way that they lead to terms resembling a 
singe field (a ^ — 1). It is important to note that even though the cancellation of higher order terms was shown only 
at linear level, the fact that it leads to a lower number of active fields in the Lagrangian means that this in fact is a 
non-perturbative effect. 

This means that apart from the case a) all the other cases lead to an effectively smaller number of active fields 
and therefore are trivial, in the sense that there exists a field redefinition that brings the Lagrangian into a function 
of a smaller number of fields. This shows that if dark energy is composed by N fields that are dynamically active 
and coupled to gravity, their presence can in principle be detected in the observables Y and rj. Needless to say, the 
effective detectability of multiple dark energy fields depends on many other factors and is beyond the scope of this 
paper. 


VI. CONCLUSIONS 


We investigated multi-field scalar-tensor models of the Brans-Dicke type, in presence of a generic source character¬ 
ized by an equation of state and a sound speed. Our aim was to find the modified gravity parameters that affect the 
Poisson equation and the anisotropic stress in the quasi-static limit. The interest in this approach lies in the fact that 
future large scale observations will probe the /c-structure of modified gravity [l^ and possibly detect whether dark 
energy is composed by one or several distinct components. 

We have shown that multi-field Brans-Dicke models with N scalar fields are described by the anisotropic stress r] 
and the effective gravitational constant Y as 


„ . Pn\k) ,, , Pi^\k) 

^ (o') ’ ^ I'Q'l 

P^"\k) Pi^\k) 


( 80 ) 
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where P^\p!?\Pn'^ are even polynomials of order n = 27V in fc, with coefficients that depend on the evolution of 
background quantities and on the matter sound speed Cg- We also find that there are no realistic non-trivial cases 
in which the higher-order k terms cancel out. Each term in the three polynomials introduces a new characteristic 
time-dependent scale kn, with n = 2,4,27V, for total of 37V scales, which can in principle be individually measured 
and employed to constrain the number of fields TV. 

An interesting question that naturally arises is how well can future surveys like e.g. the Euclid satellite [s^ [s^, 
probe the existence of more than one scalar field and constrain the values of the Brans-Dicke parameters This is 
left to future work. 
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APPENDIX 

Here we show that there always exist a single field potential V{4>) that reproduces any observed Hubble parameter 
as a function of redshift, H{z). We assume the Universe contains matter with a known generic equation of state Wm 
and a single canonical minimally coupled scalar field. The two independent Einstein equations are (we set SttG = 1) 

= Pm+ '^ + V{4>) (81) 

2H =- p^{l + Wm) (82) 

Adopting the redshift z as time coordinate, we can write from (1821) 

+ + z)H{z)'^ - Pm{z){l + Wm{z)) (83) 

where Pm{z) is the solution of 

{1 + =Zpm{z){l + Wm{z)) (84) 

From Ea. dMl) one can obtain (/)(z). Then by combining (IMll and (iMl) one obtains 

U = 3H(z)2 - (1 + z)H{z)^^ + Bl^iwmiz) - 1) (85) 

which gives V(z). Finally, by inverting (pi^z), it is possible to reconstruct V{4)) for any observed H{z). There is no 
guarantee however that the formal solution so obtained is stable or unique or free of singularities. 
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